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We study entropy production in the early stage of high-energy heavy-ion coUisions due to shear 
viscosity. We employ the second-order theory of Israel-Stewart with two different stress relaxation 
times, as appropriate for strong coupling or for a Boltzmann gas, respectively, and compare the 
hydrodynamic evolution. Based on present knowledge of initial particle production, we argue that 
entropy production is tightly constrained. We derive new limits on the shear viscosity to entropy 
density ratio t]/s, independent from elliptic flow effects, and determine the corresponding Reynolds 
number. Furthermore, we show that for a given entropy production bound, that the initial time ro 
for hydrodynamics is correlated to the viscosity. The conjectured lower bound for r]/s provides a 
lower limit for tq. 
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I. INTRODUCTION 



Experiments with colliding beams of gold ions at the 
Relativistic Heavy-Ion Collider (RHIC) have confirmed 
that dense QCD matter exhibits hydrodynamic flow ef- 
fects [ij. Their magnitude matches approximately pre- 
dictions based on ideal Euler (inviscid) hydrodynam- 
ics""^ 0. More precisely, the transverse momentum and 
centrality dependence of the azimuthally asymmetric 
flow, V2, requires a shear viscosity to entropy density 
ratio as low as 77/s < 0.2 0, i, i, 0]; this is much 
lower than perturbative extrapolations to temperatures 
T ~ 200 MeV Q. However, it is comparable to recent 
results for SU(3) pure gauge theory from the lattice 0, 
and to the conjectured lower bound for strongly coupled 
systems 77/s > l/(47r) .lOj . Similar constraints on 77/s 
have been derived from transverse momentum correla- 
tions [TT| and from energy loss and flow of heavy quarks 
at RHIC [12. 

The purpose of this paper is to obtain an independent 
upper bound on 77/s by analyzing entropy production 
in the early stages of the hydrodynamic evolution (the 
plasma phase), where the expansion rate and hence the 
entropy production rate is largest. Entropy production 
in heavy-ion collisions due to viscous effects has been 
studied before [3, [13 • The new idea pursued here is 
that recent progress in our understanding of gluon pro- 
duction in the initial state constrains the amount of ad- 
ditional entropy produced via "final-state" interactions. 



^ Numerical solutions of Euler hydrodynamics on finite grids al- 
ways involve some amount of numerical viscosity for stability. 
Reliable algorithms such as flux-corrected transport keep this 
numerical viscosity and the associated entropy production at a 
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and hence the viscosity and the thermahzation time. 
The second-order formalism for viscous hydrodynamics 
of Israel and Stewart [T^, and its application to one- 
dimensional boost-invariant Bjorken expansion [l^ . are 
briefly reviewed in section |TT1 

The initial condition for hydrodynamics, in particular 
the initial parton or entropy density in the central ra- 
pidity region, plays a crucial role. If it is close to the 
measured final-state multiplicity, this provides a strin- 
gent bound on viscous effects. The initial parton multi- 
plicity in heavy-ion collisions can, of course, not be mea- 
sured directly. Our analysis therefore necessarily relies 
on a calculation of the initial conditions (presented in 
section HIT)) . Specifically, we employ here a /cj^-factorized 
form of the "Color Glass Condensate" (CGC) approach 
which includes perturbative gluon saturation at small 
light-cone momentum fractions x [l7'|. However, differ- 
ent approaches for initial particle production, such as the 
HIJING model which relies on coUinear factorization sup- 
plemented with an additional model for the soft regime, 
also predicts multiplicities close to experiment 18] . The 
same is true when the heavy-ion collision is modeled as 
a collision of two classical Yang-Mills fields [l^- It is 
important to test these models for small systems, such 
as peripheral A -\- A (or even p + p) collisions, in order 
to constrain the entropy increase via final-state effects 
(thermahzation and viscosity). 

Section IIVI contains our main results. We show how 
the entropy production bound correlates 77/s to the ini- 
tial time for hydrodynamic evolution, tq. The entropy 
production rate grows with the expansion rate (i.e., how 
rapidly flow lines diverge from each other), and the to- 
tal amount of produced entropy is therefore rather sen- 
sitive to the early stages of the expansion. The bound 
on the viscosity depends also on the initial condition for 
the stress, which in the second-order theory is an inde- 
pendent variable and is not flxed by the viscosity and 
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the shear (unless the stress relaxation time is extremely 
short, as predicted recently from the AdS/CFT corre- 
spondence at strong coupling poj). 

In a recent paper, Lublinsky and Shuryak point out 
that if the initial time tq is assumed to be very small, that 
a resummation of the viscous corrections to all orders in 
gradients of the velocity field is required [2l!|. Here, we 
explore only the regime where tq is several times larger 
than the sound attenuation length F^, and so the stan- 
dard approach to viscous hydrodynamics should apply. 
Quantitatively, we find that an entropy bound of ~ 10% 
restricts ry/s to be at most a few times the lower bound 
iji/s = l/(47r)) conjectured from the correspondence p^ . 
On the other hand, somewhat surprisingly, we find that 
even rj/s = l/(47r) is large enough to give noticeable en- 
tropy production for thermalization times Tq well below 
1 fm/c (but still larger than Fs). Present constraints 
from initial and final multiplicities are not easily recon- 
ciled with such extremely short initial times [22i] . 

We restrict ourselves here to 1+ID Bjorken expansion. 
Given its (numerical) simplicity and the fact that the en- 
tropy production rate is largest at early times, this should 
provide a reasonable starting point. Estimates for the 
initial time tq, for the parton density and the stress at 
To, and for the viscosity to entropy density ratio rj/s are 
in fact most welcome for large-scale numerical studies 
of relativistic (second-order) dissipative fluid dynamics. 
Without guidance on the initial conditions, the hydrody- 
namic theory can at best provide qualitative results for 
heavy-ion collisions. 

We neglect any other possible source of entropy but 
shear viscosity at early times^. Even within this simpli- 
fied setting, there could be additional entropy production 
due to a viscous "hadronic corona" surrounding the fire- 
ball ^23*1, which we do not account for. Clearly, any ad- 
ditional contribution would further tighten the (upper) 
bound on ij/s and the (lower) bound on tq. We also 
assume that rj/s is constant. This does not hold over 
a very broad range of temperature [8] but should be a 
reasonable first approximation for T ~ 200-400 MeV. 

We employ natural units throughout the paper: h = 
c = ks = ^■ 



A single-component fluid is generally characterized by 
a conserved current (possibly more), iV^, the energy- 
momentum tensor T^'^ and the entropy current S^. The 
conserved quantities satisfy continuity equations. 



df.N" = , 



a^T'"' = . 



(1) 



In addition, the divergence of the entropy current has to 
be positive by the second law of thermodynamics. 



9^5'' >0. 



(2) 



For a perfect fluid, a well-defined initial-value problem 
requires the knowledge of T'^" and of iV^ on a space-like 
surface in 3+ ID Minkowski space-time. This is equiv- 
alent to specifying the initial flow field u^, the proper 
charge density n = u^N^ ^ and the proper energy density 
e = u^u^T^^; the pressure is determined via an algebraic 
relation to e and n, the equation of state (EoS). 

In dissipative fluids, irreversible viscous and heat con- 
duction processes occur. These quantities can be ex- 
pressed explicitly if the charge and entropy currents 
and the energy-momentum tensor are decomposed (pro- 
jected) into their components parallel and perpendicular 
to the flow of matter [3l|; the latter describe the dis- 
sipative currents. The transverse projector is given by 
^ _ yMy''^ with g^''' = diag(l, -1,-1,-1) the 
metric of flat space-time. In the following, we focus on lo- 
cally charge-neutral systems where all conserved currents 
vanish identically. 

The energy-momentum tensor can be decomposed in 
the following way: 

T^"' = eu^'u" -{p + H) A^'' + W^'u" + Wu^" -)- tt'^^ . (3) 

Here, = + hV^' = it^T'^^A^ is the energy 

flow, with h = {e -\- p)/n the enthalpy per particle, 
and is the heat flow; we shall define the local rest- 
frame via W^^ — (the "Landau frame"). Further- 
more, H denotes the bulk pressure such that p -I- H = 
— 'IA^^T'^^, while the symmetric and traceless part of 
the energy- momentum tensor defines the stress tensor, 



The entropy current is decomposed as 



II. DISSIPATIVE FLUID DYNAMICS 

A. Second-order formalism 

In this section we briefly review some general expres- 
sions for viscous hydrodynamics which will be useful in 
the followin g . More extensive discussions are given in 
refs. [11 m, IE HE [13, M, M, [la] , for example. 



^ This contribution is expected to vanish once transverse expansion 
is fully developed, see sect. Ill Bl 



(4) 



In the standard first order theory due to Eckart [Sll and 
Landau and Lifshitz 32] , only linear corrections are taken 
into account, i.e., = /T. On the other hand, the 
second order theory of relativistic dissipative fluid dy- 
namics includes terms to second order in the irreversible 
flows and in the stress tensor [ISl: 



5-^ = 



1^ - (/3otf - Piq.q" + /327r.„7r''") ^ 



aoHq'^ 



T 



(5) 



where the coefficients /3o, /32 and ao, oti represent 
thermodynamic integrals which (near equilibrium) are 
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related to the relaxation times of the dissipative correc- 
tions. Furthermore, from one can find linear rela- 
tionships between the thermodynamic forces and fluxes, 
leading to the transport equations describing the evolu- 
tion of dissipative flows [15| . 

In what follows, we will focus on shear effects and ne- 
glect heat flow and bulk viscosity, hence ^ simplifies 
to 

T'"' = eu'^u" - pAf"" + TT^"" . (6) 

The stress tensor satisfies a relaxation equation, 

T^u^^x^^^"' + vr^'^ ^ 27^0^'" , (7) 

where rj denotes the shear viscosity, and the shear tensor 
a*^^ is a purely "geometrical" quantity, determined by 
the flow field: 



(8) 



with ~ ISy'^dij. The relaxation time determines 
how rapidly the stress tensor tt'^" relaxes to the shear 
tensor a^^ \ in particular, in the limit — > 0, 



(9) 



satisfy the same algebraic relation as in the first-order 
theory. The limit Ttt ^ is formal, however, since the 
deviation of the stress tt'^" from 2r\(j^^ at any given time, 
as obtained by solving eq. ([7]) , depends also on its initial 
value. If ([9]) is approximately valid at the initial time 
then the first-order theory may provide a reasonable ap- 
proximation for the entire evolution (see below). 

By analyzing the correlation functions of the stress 
that lead to the definitions (|7l9p , respectively, of the shear 
viscosity, Koide argues that in the second-order theory of 
Israel and Stewart r\ may represent a different quantity 
than in the first-order approach [s^. Nevertheless, here 
we assume that the conjectured lower bound for "qjs ap- 
plies even to the causal (second-order) approach. 



B. Dissipative Bjorken scaling fluid dynamics 

In this section we recall the 1-l-ID Bjorken scaling so- 
lution [Tg'I in 3+lD space-time including stress [l3|- By 
assumption, the fiuid in the central region of a heavy-ion 
collision expands along the longitudinal z-direction only, 
with a fiow velocity v equal to z/t. This is appropriate for 
times less than the transverse diameter R of the collision 
zone divided by the speed of sound Cs = \J dp/de (possi- 
bly longer for very viscous fiuids). After that transverse 
expansion is fully developed and we expect that entropy 
production due to shear decreases. In fact, it is straight- 
forward to check that for three-dimensional scaling flow'^ 



^ Replace z —* |r| in the definition of proper time r and space-time 
rapidity f} below, and take u = r/r. 



Uf^d^a'^'^ = 0; hence, within the first-order theory at 
least, the shear viscosity does not enter the evolution 
equation of the energy density anymore. 

Formulations of the Israel-Stewart second-order theory 
for Bjorken plus transverse expansion have been pub- 
lished [H, HE m, m m, ill but require large-scale 
numerical computations. A relatively straightforward 
1+lD analysis is warranted as a first step to provide an 
estimate for entropy production. 

It is convenient to transform from (i, z) to new (r, fj) 
coordinates, where r = \Jt^ — denotes proper time 
and f]=\ log((i -I- z)/(t — z)) is the space-time rapidity; 
for the Bjorken model, it is equal to the rapidity of the 
flow, 77fi = ^log((l -I- v)l{\ — v)). In other words, the 
four- velocity of the fiuid is = (5^*^ -h ^^'^) x^'It. 

The longitudinal projection of the continuity equation 
for the stress-energy tensor then yields 



de 
rf7 



V 



T 



n 

r 



0. 



(10) 



Here, e is the energy density of the fiuid in the local rest- 
frame, while p denotes the pressure. These quantities 
are related through the equation of state (EoS) . We focus 
here on entropy production during the early stages of the 
evolution where the temperature is larger than the QCD 
cross-over temperature T^— 170 MeV, and so assume a 
simple ideal-gas EoS, p = e/3. 

In what follows, we will neglect the bulk pressure 11 
which would otherwise tend to increase entropy produc- 
tion further. Well above Tc, this contribution is expected 
to be much smaller than that due to shear t34|. In 
the transition region, the bulk viscosity could be signifi- 
cant [1^. Also, for 1-l-lD expansion considered here, the 
stress 4> = tt'^" — tt^^ acts in the same way as the bulk 
pressure 11: only the combination $ — 11 appears in (jlOJ. 

The time evolution of the stress is determined by (14| 





$ 


'1 


T 


d 


dr 


^2 


T 




dr 



= 



(11) 



sets the time scale for relaxation to the first- order 
theory where 4>ist-o — 477/3r (not to the ideal-fiuid limit 
<f> = 0). It is related to 77 and (32 via r^r = 2ry/32- For a clas- 
sical Boltzmann gas of massless particles P2 = 3/(4p) = 
3/(Ts) and so 
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(12) 



At infinite coupling, from the AdS/CFT correspon- 
dence r]/s = l/(47r) [10] and = (1 - log2)/(67rT) 
For large but finite coupling, we assume that ril32 and 
thus are proportional to ry/s, i.e. that /?2 = (37')/(Ts) 
with r (1 - log2)/9; then 
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(13) 



Note that the numerical prefactor r is about 30 times 
smaller than for a Boltzmann gas, implying much faster 
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relaxation of the dissipative fluxes to the first-order the- 
ory. 

One can also define a Reynolds number via the ratio 
of non-dissipative to dissipative quantities (36j . i? = (e + 
p)/<i>. Eq. ([To]) can then be written as 



lower bound on tq: 



de 



d log T 



(i?-i-l)(e + p). 



(14) 



where we have neglected the bulk pressure. For stability, 
the effective enthalpy {e+p){l — 1/i?) should be positive, 
i.e. R > 1. The energy density then decreases monotoni- 
cally with time. 

The equations of second-order dissipative fluid dynam- 
ics, (Ho]) or dH]) and (HI]) together with ^ or ^ and 
f32 = TTr/{2ri) form a closed set of equations for a fluid 
with vanishing currents, if augmented by an EoS. Fur- 
thermore, the initial energy density cq = e(T'o) f^nd the 
initial shear $o = ^{tq) have to be given. In the second- 
order theory one has to specify the initial condition for 
the viscous stress $o independently from the initial en- 
ergy or particle density. We are presently unable to com- 
pute $0. Below, we shall therefore present results for 
various values of $0 • 

Alternatively, a physically motivated initial value <i>g 
for the stress can be obtained from the condition that 
dR/dr = at T = Tq. This is the "tipping point" between 
a system that is already approaching perfect fluidity at 
To (dR^^/dr < if $0 > 'i'o) ^^'^ oi^^ that is unable to 
compete with the expansion and is in fact departing from 
equilibrium {dR~^/dT > if $0 < $0) for at least some 
time after tq. 

For an EoS with constant speed of sound, say p = e/3, 
the condition that R = is equivalent to e/e = $/$. 
Eqs. (dn]) and ^ then yield 



n 

eo 
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(Ist-order theory) 



(15) 

(16) 
(17) 



The second line applies in the limit of short relaxation 
time; since t^^ is proportional to r, this is always satisfied 
in the limit r ^ 0. For typical initial conditions relevant 
for heavy-ion collisions, it is a reasonable approximation 
even in the Boltzmann limit (r = 1). In (|17p we have 
indicated that (I16|) is in fact nothing but the stress in the 
first-order theory (divided by the initial energy density). 
While it is clear that $ relaxes to <i>ist-o over time-scales 
on the order of t^^ , eq. (|17p is actually a statement about 
the initial value of $: in a fluid with reasonably short 
relaxation time and stationary initial Reynolds number, 
R-iTo) = 0, even the initial value of the stress is given by 
the first-order approach. 

The condition i?(To) > 1 for applicability of hydrody- 
namics together with eq. (|16p then provides the following 
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To > -- — = Fs To , 
3To s 



(18) 



where F^ denotes the sound attenuation length; within 
the first-order approach, R = t/Ts- The factor of rj/s 
on the right-hand-side illustrates the extended range of 
applicability of hydrodynamics as compared to a Boltz- 
mann equation: a classical Boltzmann description re- 
quires that the thermal de-Broglie wave length, ~ 1/T, is 
smaller than the (longitudinal) size of the system, t. For 
very small viscosity, though, hydrodynamics is applicable 
(since i? ^ 1) even when F^ ^ t ^ 1/T . In this point we 
differ somewhat from Lublinsky and Shuryak |2l| , who 
argue that the theory needs to be resummed to all or- 
ders in the gradients of the velocity field already when 
T ~ 1/T. From our argument above, this should be nec- 
essary only when To ~ Fs(to), which is much smaller than 
1/To if r?/s < 1. 

The purpose of this paper is to motivate, however, that 
a much stronger constraint than To > Fs(to) may anyhow 
result from a bound on entropy production (which follows 
from the centrality dependence of the multiplicity), cf. 
section IIVI 

In the Bjorken model, the entropy per unit rapidity 
and transverse area at time t is given by 



1 dS{T 
d-q 



ts{t) , 



(19) 



where A± is the transverse area while s = S^u^ denotes 
the longitudinal projection of the entropy current. Ne- 
glecting heat flow {q^ = 0) and bulk pressure (Ft = 0) 
one obtains from ([5]): 



4rs 



(20) 



s can be determined, for any t > to, from the solution 
of eqs. (|10llip . Note that the second term in (PO)) is of 
order ($/e)^. For nearly perfect fluids with rj/s ^ 1 and 
R{tq) = it is rather small. 



III. THE CGC INITIAL CONDITION 

Before we can present solutions of the hydrodynamic 
equations, we need to determine suitable initial condi- 
tions. To date, the most successful description of the 
centrality dependence of the multiplicity is provided by 
the Kharzeev-Levin-Nardi (KLN) fcj^-factorization ap- 
proach [17] . The KLN ansatz for the unintegrated gluon 
distribution functions (uGDF) of the colliding nuclei in- 
corporates perturbative gluon saturation at high energies 
and determines the p^-integrated multiplicity from weak- 
coupling QCD without additional models for soft particle 
production. 
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Specifically, the number of gluons that are released 
from the wavefunctions of the colliding nuclei is given 
by 



(Pr^dy 



N, 



Pi 



X(j)A{xi,{p± + k±)^ /4:;ri_) 
X(j)B{x2, (p±-fcj_)^/4;r^) 



(21) 



where TVc = 3 is the number of colors, and p±, y are the 
transverse momentum and the rapidity of the produced 
gluons, respectively. xi^2 = P± exp(±i; ) /^/s n n denote 
the light-cone momentum fractions of the colliding gluon 
ladders, and ^/snn — 200 GeV is the collision energy. 

The normalization factor A/" can be fixed from periph- 
eral collisions, where final-state interactions should be 
suppressed. (Ideally, the normalization could be fixed 
bom p+p collisions; however, this is possible only at suf- 
ficiently high energies, when the proton saturation scale 
is at least a few times Aqcd-) A/" also absorbs NLO cor- 
rections; when we compare to measured multiplicities of 
charged hadrons, it includes as well a factor for the av- 
erage charged hadron multiplicity per gluon, and a Jaco- 
bian for the conversion from rapidity to pseudo-rapidity. 

The uGDFs are written as follows [l^,!!^,!!!: 



(x, k I 



oisiQV) niax(Q2^A:2 



■P{r^){\-x)' 



(22) 

P(r^) denotes t he p robability of finding at least one nu- 
cleon at r±_ [s^, H^] . This factor arises because config- 
urations without a nucleon at r±^ do not contribute to 
particle production. Note that the perturbative ^ \/k\ 
growth of the gluon density towards small transverse 
momentum saturates at k± — Qg. Therefore, the p±- 
integrated gluon multiplicity obtained from (PTjl is finite. 

We should emphasize that the ansatz is too sim- 
ple for an accurate description of high-pj^ particle pro- 
duction. For example, it does not incorporate the so- 
called "extended geometric scaling" regime above Qs, 
which plays an important role in our understanding of 
the evolution of high-p^ spectra from mid- to forward 
rapidity in d-|-Au collisions [s^. However, high-pj^ par- 
ticles contribute little to the total multiplicity, and more 
sophisticated models for the uGDF do not change the 
centrality dependence of dN/ dy significantly |32| ■ 

Qsix, r±) denotes the saturation momentum at a given 
momentum fraction x and transverse coordinate r±. It 
is parameterized as 37, 38] 



g2(x,r-i) = 2GeV- 



, / r(rJ/P(r^) \ 



1.53 



0.01 



(23) 



The ~ l/x^ growth at small x is expected from BFKL 
evolution and has been verified both in deep inelastic 
scattering at HERA and in high-p i particle produc- 
tion from d + Au collisions at RHIC [33|; the growth 
speed is approximately A ~ 0.28. Note that the satura- 
tion momentum, as defined in (12311 . is "universal" in that 



it doesn't depend on the thickness of the collision partner 
at r±_ [4l|. 

The centrality dependence of Qs is determined by the 
thickness function T(r^), which is simply the density 
distribution of a nucleus, integrated over the longitudi- 
nal coordinate z. Note that the standard Woods-Saxon 
density distribution is averaged over all nucleon configu- 
rations, including those without any nucleon at r j. For 
this reason, a factor of 1/P(r^) arises in Ql [33,[38|. It 
prevents Qa from dropping to arbitrarily small values at 
the surface of a nucleus (since at least one nucleon must 
be present at r_\_ or else no gluon is produced at that 
point). The fact that Qs is bound from below prevents 
infrared sensitive contributions from the surface of the 
nucleus and also makes the uGDF p2|) less dependent 
on "freezing" of the one-loop running coupling. 
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FIG. 1: (Color online) Centrality dependence of the charged 
particle multiplicity at midrapidity from the fcx-factorization 
approach with perturbative gluon saturation at small-x, for 
Cu-f Cu and Au+Au collisions at full RHIC energy, ^/snn = 
200 GeV. PHOBOS data from ref. ^\ the errors are system- 
atic, not statistical. 

Fig. [1] shows the centrality dependence of the mul- 
tiplicity, as obtained from eq. \21\ via an integration 
over the transverse plane. It was noted in 37] that 
the multiplicity in the most central collisions is signifi- 
cantly closer to the data than the original KLN predic- 
tion if the integration over r j_ is performed explicitly, 
rather than employing a mean-field like approximation, 
QH^i^) (<3s)(^)- All sven better description of the 
data can be obtained when event-by-event fluctuations of 
the positions of the nucleons are taken into account [ssj ; 
they lead to a slightly steeper centrality dependence of 
the multiplicity per participant for very peripheral colli- 
sions or small nuclei. However, we focus here on central 
Au-|-Au collisions and hence we neglect this effect. 

It is clear from the figure that the above CGC-fc_L- 
factorization approach does not leave a lot of room for an 
additional centrality- dependent contribution to the parti- 
cle multiplicity. (Centrality independent gluon multipli- 
cation processes have been absorbed into M .) In fact. 
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within the bottom-up thermahzation scenario one does 
expect, parametrically, that gluon splittings increase the 
multiphcity by a factor ^ Xjc?!^ [4^ before the system 
thermahzes at tq and the hydrodynamic evolution be- 
gins. If the scale for running of the coupling is set by 
Qs, this would lead to an increase of the multiplicity 
for the most central Au-|-Au collisions by roughly 20%. 
However, such a contribution does not seem to be vis- 
ible in the RHIC data, perhaps because the bottom-up 
scenario, which considered asymptotic energies, does not 
apply quantitatively at RHIC energy. It is also think- 
able that the model (|22l23p overpredicts the growth of 
the particle multiplicity per participant with centrality 
somewhat. 

It is noteworthy that from the most peripheral Cu+Cu 
to the most central Au-|-Au bin, that [dN / drf) / Npa,^t 

1 /3 

grows by only ~ 50% while N^'^^^ increases by a factor of 
2.6. Clearly, any particle production model that includes 
a substantial contribution from perturbative QCD pro- 
cesses will cover most of the growth. This implies that 
rather little entropy production appears to occur after 
the initial radiation field decoheres. If so, this allows us 
to correlate the thermahzation time tq and the viscos- 
ity to entropy density ratio -q/s. We shall assume that 
about 10% entropy production may be allowed for central 
Au-f Au collisions. 

The density of gluons at = 1/Qs is given by 
dN/d^r±dy from eq. (^1]) . divided by Tg. For a central 
collision of Au nuclei at full RHIC energy, the average 
Qs — 1.4 GeV at midrapidity; hence Tg ~ 0.14 fm/c. The 
parton density at this time is approximately ~ 40 fm~'^. 
If their number is effectively conserved^ until thermahza- 
tion at To then 
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FIG. 2: (Color online) The parton number and energy den- 
sities (averaged over the transverse plane) for 6 = Au+Au 
collisions at y^sjvjv ~ 200 GeV as functions of the thermal- 
ization time tq. 



IV. RESULTS 



A. Evolution of the entropy and of the Reynolds 
number 



We begin by illustrating entropy production due to dis- 
sipative effects. Given an initial time tq for the hydrody- 
namic evolution, we determine AS = dStin/ dy — dS-mi/ dy 
for T > Tq. This quantity increases rather rapidly at first, 
since the expansion rate H = d^u^ = l/r is biggest at 
small T. We chose Tfin = 5 fm/c to be on the order of 
the radius of the collision zone and fix the final value for 
AS/Sini to equal 10%. Having fixed the initial and final 
entropy as well as the initial time then determines t]/s. 



"(to) 



To 



■ n{Ts) . 



(24) 



The initial energy density e(To) can now be obtained from 
the density via standard thermodynamic relations. We 
assume that the energy density corresponds to 16 gluons 
and 3 massless quark flavors in chemical equilibrium. 



2^4 _ 43C|3)^3 ^25) 



e(r) = ^.^r^ , n{T) 

30 TT 



Fig. [2] shows the parton number and energy densities at 

To. 



* We repeat that centrality-independent gluon muhipHcation pro- 
cesses and the contribution from quarks are already accounted 
for via the normaUzation factor Af. 



C/0 
<1 



0.1 



0.08 



.S 0.06 



0.04 
0.02 




rj/s = 0.035 
rj/s = 0.066 
ri/s = 0.105 



1 2 3 

X (fm/c) 



FIG. 3: (Color online) Entropy production within second- 
order dissipative hydrodynamics as a function of (proper) 
time; the initial value of the stress is $(ro) — $oi cf. 
eq. (|15p . Curves for three different initial times, to = 0.5, 
1, 1.5 fm/c are shown, and for each curve 77/s is chosen such 
that AS/Sini = 10% at tr^ = 5 fm/c. 
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The result of the calculation is shown in Fig. [3l As ex- 
pected, if the hydrodynamic expansion starts later (larger 
To) then less entropy is produced for a given value of 
rj/ s] conversely, for fixed entropy increase, larger values 
of 77/5 are possible. This is due to two reasons: the to- 
tal time interval for one-dimensional hydrodynamic ex- 
pansion as well as the entropy production rate decrease. 
In fact, the figure shows that for very small initial time 
the AS'/S'ini = 10% bound can not be satisfied with 
ri/s > l/(47r) ~ 0.08. 



0.15 




0.05 



1.5 2 2.5 3 3.5 4 
X (fm/c) 

FIG. 4: (Color online) Time evolution of the inverse Reynolds 
number (using the Boltzmann relaxation time) for different 
initial values of the viscous stress, $0, and of the viscosity 
to entropy density ratio ri/s. The initial time is ro = 1 fm/c 
for all curves. The short-dashed line corresponds to the ini- 
tial condition R{to) = 0, cf. eq. (|15p . The long-dashed line 
corresponds to the first-order theory. 



In Fig. 3] we show the behavior of the inverse Reynolds 
number for different initial values of the stress. Again, for 
each curve 77/s is fixed such that AS/ Sini = 10% at Tfi„ = 
5 fm/c. As already indicated above, if <I>o < $0 defined 



in eq. (jlSp , the fluid can not compete with the expansion 
and departs from equilibrium. On the other hand, if $0 > 
<I>Q , there is already a rapid approach towards the perfect- 
fluid limit at tq. In either case, the interpretation of tq 
as the earliest possible starting time for hydrodynamic 
evolution does not appear sensible. The initial condition 
corresponding to R{to) = in turn corresponds to the 
situation where the fluid has just reached the ability to 
approach equilibrium. It is clear from the figure that 
the evolution is close to that predicted by the first-order 
theory. 

Fig. [S] shows the Reynolds number for our ansatz 
for the relaxation time at strong coupling, which essen- 
tially follows the behavior given by the first-order theory: 
after a time ~ has elapsed, the behavior of R is nearly 
independent of the initial value of $. The initial condi- 
tion <I>o = $0 again leads to the most natural behavior 
of R without a very rapid initial evolution. 



0.2 



0.15 



0.1 



0.05 



<I>o = 0.0 eg; Tj/s = 0.082 
<I>*0 = 0.1 18 r\ls = 0.08 



0.2 eg; Tj/s = 0.079 






0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
X (fm/c) 

FIG. 5: (Color online) Same as Fig. |3]but for the generalized 
AdS/CFT relaxation time (fT3)) . 



B. 77/s versus tq 

From the previous results it is evident that fixing the 
amount of produced entropy, AS'/5ini, correlates 77/s 
with Tq. In this section we show the upper limit of 77/s 
as a function of tq. 



0.25 



0.2 



X 0.15 

eg 

B 



0.1 



0.05 



<l>Q = O.Oeo 

<I>o = 0.1eo --- 
■ <I>o = 0-2eo 


1 

J 1 
/ 1 

/ ' 
/ / 


1/(471) / 


/ / 
/ / 
/ / 
/ / 
/ ' 
/ / 
/ 

/ 
r 





0.5 



1 1.5 

Xq (fm/c) 



2.5 



FIG. 6: (Color online) Bound on 77/s as a function of tq, for 
AS/Sini = 10% entropy production, for a Boltzmann gas. 

We begin with the Boltzmann gas with fixed $0 (inde- 
pendent of To) in Fig. [6] One observes that the maximal 
viscosity depends rather strongly on the initial value of 
the stress. For any given $0i (?7/s)max first grows approx- 
imately linearly with tq. For large initial time, however, 
the expansion and entropy production rates drop so much 
that the bound on viscosity eventually disappears. Fur- 
thermore, it is interesting to observe that the conjectured 
lower bound r7/s = l/(47r) excludes too rapid thermal- 
ization: even if the fiuid is initially perfectly equilibrated 
('I'o = 0), a thermalization time well below ~ 1 fm/c is 
possible only if either r7/s < l/(47r) or AS'/5'ini > 10%. 
With 10% corrections to perfect fluidity at tq, shown 
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by the long-dashed hne in Fig. [6l the minimal tq com- 
patible with both 77/s > l/(47r) and AS'/5i„i = 10% is 
about 1.2 fm/c. If rj/s ~ 0.1 — 0.2, as deduced from the 
centrality dependence of elliptic flow at RHIC then 
To — 1.5 fm/c. 



0.25 



0.15 




0.05 



0.5 1 1.5 

Xq (fm/c) 



2.5 



FIG. 7: (Color online) Same as Fig.|6]but for strong coupling, 
Tn from (fT5)) . 

In Fig. [7] we perform a similar analysis for our 
ansatz for the strong-coupling case. Due to the 
much smaller relaxation time of the viscous stress, we ob- 
serve that the viscosity bound is now rather insensitive 
to the magnitude of the initial correction to equilibrium. 
We obtain a lower bound on the thermalization time of 
To — 1 fm/c for the minimal viscosity, increasing to about 
1.2 - 2.2 fm/c if r^/s ~ 0.1 - 0.2. 



0.25 



0.2 - 



AS/Sini= 5% 
AS/Sini = 1 0% 
AS/S;„; = 1 5% 




0.5 



1 1.5 

Xq (fm/c) 



2.5 



FIG. 8: (Color online) Same as Fig. |6] but for $0 = 'I'o and 
different bounds on entropy production. 

In Fig. [8] we return to the Boltzmann gas with initial 
value for the stress as given in eq. (jl5p . corresponding 
to R{tq) = 0. Comparing to Fig. [6l we observe that the 
viscosity bound is affected mostly for large To: with this 
initial condition, a high viscosity rj/s ^ 1 is excluded even 
if the initial time is as big as 2 fm/c. The reason why the 



upper bound on the viscosity does not disappear at large 
Tq for this initial stress is that ^q/cq grows with 77/s, cf. 
eq. (fT5|). A lot of entropy would then be produced, even 
for large To. 
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FIG. 9; (Color online) Same as Fig. [8] but for the strong- 
coupling limit, i.e. for Ttt from eq. (I13p . 



Fig. [8] also gives an impression of the sensitivity to the 
entropy production bound. For 77/s = 0.15, for example. 
To decreases from ~ 1.8 fm/c, if the entropy is allowed to 
increase by 10%, to ~ 1.5 fm/c if the bound is relaxed to 
AS'/5'i„i = 15%. 

We performed similar calculations for the strong- 
coupling limit as shown in Fig. [9l The curves are rather 
close to those for a Boltzmann gas from Fig. [51 which is 
expected. With this initial condition, i.e. $0 = ^oi 
hydrodynamic evolution is close to the first-order theory 
for both cases. Entropy production is sensitive only to tq 
and 7;/s but is nearly independent of the stress relaxation 
time T-jr. 



0.8 



00 



0.6 



.5 0.4 



00 



0.2 




Boltzmann; r|/s = 1/(4;t) 
AdS/CFT; ti/s = 1/(47t) 

Boltzmann; t|/s = 1/(2it) 
AdS/CFT; ti/s = 1/(27t) 



0.5 
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FIG. 10: (Color online) Initial over final entropy as a func- 
tion of To for two different values of 77/s, and two different 
relaxation times corresponding to weak and strong coupling. 

Fig. [TUl shows the ratio of the final to the initial entropy 
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as a function of tq for two different values of rj/s, and the 
two different relaxation times discussed above in eqs. 
fT3| . Here, S'fin has been fixed to the value appropriate for 
central Au+Au collisions while S'ini is varied accordingly. 
For example, for tq — 0.6 fm/c and rj/s — l/(27r), almost 
30% entropy production occurs. This would account for 
the entire growth of {dN /drj)/N ■ps.^t from iVpart — 60 to 
A'part — 360 observed in Fig. [TJ That is, for these pa- 
rameters the initial parton multiplicity per participant 
would have to be completely independent of centrality. 
For the same viscosity, tq = 0.3 fm/c would imply that 
nearly half of the final-state entropy was produced dur- 
ing the hydrodynamic stage, i.e. that the initial multi- 
plicity per participant should actually decrease with cen- 
trality. Such a scenario appears unlikely to us. Note 
that even for tq = 0.3 fm/c and -q/s — l/(47r), with 
T = 400 MeV one finds that Vs{jq)Itq ~ 0.17 is quite 
small. Romatschke obtained similar numbers for the ini- 
tial to final entropy ratio, albeit only for tq = 1 fm/c, 
in a computation that included cylindrically symmetric 
transverse expansion (28j . 



V. SUMMARY, DISCUSSION AND OUTLOOK 

In this paper, we have analyzed entropy production 
due to non-zero shear viscosity in central Au-I-Au col- 
lisions at RHIC. We point out that a good knowledge 
of the initial conditions, and of the final state, of course, 
can provide useful constraints for hydrodynamics of high- 
energy collisions, specifically on transport coefficients, on 
the equation of state (not discussed here, cf. [Q]), on the 
initial/thermalization time and so on. 

Our main results are as follows. Assuming that hy- 
drodynamics applies at T > Fs, then due to the rather 
restrictive bound on entropy production, it follows that 
the shear viscosity to entropy density ratio of the QCD 
matter produced at central rapidity should be small, at 
most a few times the lower bound rj/s = l/(47r) conjec- 
tured from the AdS/CFT correspondence at infinite cou- 
pling. This represents a consistency-check with similar 
numbers {rj/s < 0.2) extracted from azimuthally asym- 
metric elliptic flow 0, H, [1, 0] ■ We have neglected sev- 
eral other possible sources of entropy production, such as 
bulk viscosity near the transition region [s^ or hadronic 
corona effects [23| : such additional contributions might 
tighten the constraints even further. 

Furthermore, the entropy production bound correlates 
the maximal allowed viscosity to the initial time tq for 
hydrodynamic evolution. This is due to the fact that the 
expansion rate is equal to the inverse of the expansion 
time, which makes entropy production from viscous ef- 
fects rather sensitive to the value of tq. We have found 
that for A5'/S'ini — 10%, that the initial time for hy- 
drodynamics should be around 1 fm/c, possibly a little 
larger. Significantly smaller thermalization times would 
either require rj/s < 0.15 — 0.2 (or even smaller than 
l/(47r)). Alternatively, they would require a particle pro- 



duction mechanism that yields significantly lower initial 
multiplicities than the KLN-CGC approach. Given the 
very good description of the centrality dependence of the 
multiplicity, however, to us it appears reasonable to as- 
sume that this approach provides an adequate initial con- 
dition in that the initial parton multiplicity per partici- 
pant increases with centrality. 

A significant problem with viscous hydrodynamics, 
in particular with the second-order approach of Israel- 
Stewart, is the fact that the number of initial parameters 
increases. Even within the most simple framework fol- 
lowed here (1-l-lD Bjorken expansion combined with ne- 
glect of conserved currents, of bulk viscosity, and of heat 
flow) , a unique solution requires us to specify, in addition 
to the ideal-fluid parameters, the shear viscosity and the 
initial value for the stress. The latter, in particular, is not 
a general property of near-equilibrium QCD but depends 
on the parton liberation and thermalization process. We 
have, however, introduced a physically motivated initial 
condition for the stress: if tq is defined as the earliest 
possible initial time for hydrodynamics, it is plausible 
that the initial Reynolds number should be stationary, 
R{to) = 0. Otherwise the fluid either still departs from 
equilibrium (R(tq) < 0); or is already approaching it 
{R{ro) > 0). 

For small relaxation times of the stress, the condition 
that R{to) = implies that its initial value already be 
close to that given by the first-order theory of Eckart, 
Landau and Lifshitz (the relativistic generalization of 
Navier-Stokes hydrodynamics). We therefore expect that 
in general the two approaches will provide rather similar 
results for heavy-ion collisions. One should keep in mind, 
however, that in the second-order theory the entropy cur- 
rent includes a term quadratic in the stress, which is of 
course absent from the first-order theory, and which re- 
duces entropy production slightly. 

Perhaps most importantly, with R{to) — 0, the hy- 
drodynamic evolution is largely independent of the stress 
relaxation time , and therefore similar for both a Boltz- 
mann gas at weak coupling (with low viscosity, however) 
and a strongly coupled plasma. The latter relaxes very 
rapidly to the first-order theory, regardless of the initial 
condition. The former, on the other hand, is forced by 
the initial condition to start close to relativistic Navier- 
Stokes, and the relaxation time is still sufficiently small to 
prevent a significant departure from the first-order the- 
ory. 

The initial condition R{to) = also guarantees that 
R{t) ^ 1 for all r > tq, as long as the initial time 
is not extremely short {tqTq ^ rj/s). The effective en- 
thalpy (1 — 1/R){e + p) is therefore always positive. On 
the other hand, our numerical results indicate that the 
Reynolds number does not exceed ~ 100 during the QGP 
phase. This is well below the regime where Navier-Stokes 
turbulence occurs in incompressible, non-relativistic flu- 
ids {R <; 1000). Indeed, turbulence during the hydrody- 
namic stage would probably cause large fluctuations of 
the elliptic flow V2 [Si , which are not seen [i^ . 
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A quantitative interpretation of hydrodynamic flow ef- 
fects in heavy-ion collisions at RHIC and LHC will of 
course require 2-MD and 3-l-lD solutions 0, [13, [H, S 
H^jlll, [23- The results obtained here should prove use- 
ful for constraining the initial conditions (in particular 
To and <i>o) for such large-scale numerical efforts. In par- 
ticular, as we pointed out here, the entropy production 
bound correlates tq with rj/s. In turn, we expect that 
elliptic flow will provide an anti- correlation since later 
times and larger shear viscosity should both reduce its 
magnitude. The intersection of those curves could then 
provide an estimate of the initial time for hydrodynamics 



at RHIC. 
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